Chemical bonding, elasticity, and valence force field models: a case study for Q:-Pt2Si 

and PtSi 
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We have carried out a detailed study of the chemical bonding for two room-temperature stable 
platinum silicide phases, tetragonal Q-Pt2Si and orthorhombic PtSi. An analysis of the valence 
electronic charge density reveals surprising evidence of covalent three-center bonds in both silicide 
phases, as well as two-dimensional metallic sheets in a-Pt2Si. These elements of the bonding are fur- 
ther analyzed by constructing valence force field models using the results from recent first principles 
calculations of the six (nine) independent, non-zero elastic constants of a-Pt2Si (PtSi). The resulting 
volume-, radial-, and angular-dependent force constants provide insight into the relative strength of 
various bonding elements as well as the trends observed in the elastic constants themselves. The 
valence force field analysis yields quantitative information about the nature of the chemical bonding 
which is not easily discernable from the more qualitative charge density plots. More generally, this 
study demonstrates that the detailed variations in the elastic constants of a material contain use- 
ful information about the chemical bonds which can be extracted using valence force field models. 
Inversely, these models also allow identification of specific elements of the chemical bonding with 
particular trends in the elastic constants, both within a given material and among a class of related 
materials. 

PACS 61.50.Lt, 62.20.Dc 



I. INTRODUCTION 



Deposition of metallic platinum silicide compounds 
on silicon substrates leads to the formation of rectify- 
ing junctions, with a Schottky barrier of 220-240 meV 
(for liolcaLin the case of orthorhombic PtSi on p-type 
Si (OOl).El'0 This energy matches an important atmo- 
spheric "transparency window" in the infrared region, 
making these materials well suited to infrared detector 
applications. PtSi has also been discussed as a promis- 
ing candidate to replace Ti2Si in polysilicon ialjfirconnect 
applications in sub-half-micron technologies JjU In light 
of these and other technological applications, as well as a 
general paucity of earlier treatments of the fundamental 
properties of the platinum silicides, there have been two 
recent in-depth studies of the atomic and electronic struc- 
ture of two room-temperature stable platinum silicide 
phases, orthorhombic PtSi and tetragonal Q;-Pt2Si. Beck- 
stein et al.a have carried out an extensive set of first prin- 
ciples electronic structure calculations for both materials. 
In addition to the electronic structure, they have calcu- 
lated all of the equilibrium structural parameters and 
zero-gressure elastic constants for both phases. Franco 
et al.u used a combination of photoelectron spectroscopy 
(PES), soft x-ray emission spectroscopy (SXE), and x- 
ray absorption spectroscopy (XAS) to study the detailed 
electronic structure of orthorhombic PtSi. First princi- 
ples calculations of the partial density of states (PDOS) 
were also carried out in order to aid in interpreting the 
experimental spectra. 

The present study is complimentary to these two ear- 



lier treatments and makes contact with them in a number 
of ways. The combination of the atomic and electronic 
structure gives rise to the chemical bonding of a mate- 
rial. The elastic constants and the various experimental 
spectroscopies reflect the details of this bonding but they 
do so indirectly. One of the goals here is to directly eluci- 
date the fundamental nature of the chemical bonding in 
the two silicide phases studied previously. Towards that 
end we have calculated and analyzed the valence elec- 
tronic charge density for both silicides. However, this 
analysis is only qualitative and thus we have made fur- 
ther attempts to gain a more quantitative understanding. 
The previous first principles study noted a number of in- 
teresting trends in the elastic constants, both within a 
given material and amoag the two silicides and the pure 
Pt and pure Si phases.tl In the present work we ana- 
lyze these trends in much greater detail and in a more 
quantitative fashion by constructing valence force field 
models for all four materials. The models are obtained 
by fitting the first principles elastic constants while also 
using insights gained from the charge density analysis to 
guide the particular choice of radial and angular inter- 
actions. In turn, the magnitudes of these various inter- 
actions, as obtained from the fits, provide a quantitative 
measure of the relative importance of different elements 
of the chemical bonding. In addition, the models can 
be inverted by expressing the various elastic constants 
in terms of the volume-, radial-, and angular-dependent 
interactions. We are thus able to identify the individual 
trends in the elastic constants with particular elements 
of the chemical bonding. 
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In the present work we have two overall goals. The first 
is to gain a quantitative understanding of the chemical 
bonding in tetragonal Q!-Pt2Si and orthorhombic PtSi. 
The second goal is to demonstrate, through a case study 
of these two silicides as well as pure Pt and pure Si, that 
in general terms the variations of the elastic constants 
of a material contain useful information about the chem- 
ical bonding and that valence force field models are a 
convenient means for extracting this information. More- 
over, by inverting the models and identifying the chemi- 
cal interactions responsible for the observed trends in the 
elastic constants we are thereby able to obtain a more 
intuitive understanding of the connection between chem- 
ical bonding and the mechanical properties of a material. 
Given this more general goal we have therefore described 
the construction of the models and the analysis of var- 
ious elastic constant trends in some detail. Sec. || pro- 
vides the relevant details regarding the atomic structure 
for the two platinum silicides studied here and in Sec. |l| 
we summarize the previous elastic constant calculations 
from Ref. ^. The valence electronic charge densities are 
analyzed in Sec. |^ and the valence force field models 
are presented in Sec. Our results are summarized in 



Sec. VI 



II. ATOMIC STRUCTURE 



The stable phase of puse Pt at ambient conditions is 
face-centered cubic (fcc),El while for pure Si it is cubic 
diamondJj The conventional unit cells of the two plat- 
inum silicides a-Pt2Si and PtSi are shown in Fig. 0. The 
room-temperature (T < 968 K) a-phase of PtjSi|-occurs 
in the body-centered tetragonal (bet) structureJl^'y The 
Strukturbericht desigaation is i'2f, and the space group is 
lA/mmm (No. 139). 113 The two symmetry-equivalent Pt 
atoms in the primitive cell occupy Wyckoff 4(d) sites and 
the one Si atom occupies a 2(a) site. The atom positions 
are completely determined by the space group symmetry 
but there are two independent lattice constants a and 
c. PtSi has a primitive orthorhombic structure with four 
symmetry-equivalent Pt atoms occupying Wyckoff 4(c) 
sites and four sjpaoietry-equivalent Si atoms also occu- 
pying 4(c) sites.E5~tJ The Strukturbericht designation 
for this MnP-tyjie lattice is i?31 and the space group is 
Pnma (No. 62). E3 The atom coordinates along the a- and 
c-axes are not completely specified by the space group 
symmetry and thus there are four free internal structural 
parameters upt, Wpt, usi, and v^i- The structure also has 
three independent lattice constants a, 6, and c. All of 
the relevant equilibrium structural parameters for each 
of these four materials are given in Table |, including 
both the experimental values and the self-consistent the- 
oretical values calculated from first principles in Ref. ||. 



III. ELASTIC CONSTANTS 

Since we will rely heavily on the detailed results of the 
first principles elastic constant calculations from Ref. ||, 
we summarize them here. The internal energy E of the 
crystal is expanded to second order in the elements of the 
strain tensor e^, using Voigt notation, 

V ■ 



E{V, {ej) = E{V) + vY^a.e, + - ^ c.,e.e. 



(1) 



where V is the volume of the unstrained crystal, the 
Ui are the components of the applied stress terisor, and 
the Cij are the second order elastic constantsfl Since 
the undistorted crystal was always taken to be the zero- 
pressure theoretical equilibrium structure, the applied 
stress components <Ti are all zero and so the second term 
of Eq. (|^) did not enter. All of the elastic constant calcu- 
lations were carried out using the theoretical equilibrium 
structural parameters listed in Table ^. 

Crystals with cubic space group symmetry have only 
three distinct, non-vanishing elastic constants. The the- 
oretical values of these three elastic constants for both 
pure Pt and pure Si, as obtained in Ref. ^ are listed in 
Table || together with the corresponding experimental 
values. The theoretical bulk moduli were obtained from 
the theoretical elastic constants [B^ — i(cii -|-2ci2)]. We 
note that in the case of Si the calculation of C44 required a 
relaxation of the positions of the Si atoms within the dis- 
torted unit cell. The C44 strain-induced symmetry reduc- 
tion allowed the Si atoms to relax in the [001] direction. 
Without this relaxation, C44 would have been 108.6 GPa; 
this is to be compared with the relaxed vakie of 79.9 GPa 
and the experimental value of 79.1 GPa.ll3 The require- 
ment of mechanical stability in a cubic crystal leads to 
the following restrictions on the elastic constantal3 

(cii - C12) > 0, cii > 0, C44 > 0, (cii + 2ci2) > 0. 

(2) 

These stability conditions also lead to a restriction on the 
magnitude of Bq (Ref. 0) 



C12 < -Bo < cii. 



(3) 



Tetragonal a-Pt2Si has six independent and non-zero 
elastic constants. Three of these elastic constants, cn, 
C12, and C44, correspond to strain-induced symmetry re- 
ductions for which the positions of the Pt atoms are no 
longer completely fixed by the symmetry. The strain- 
induced forces drive them into energetically more favor- 
able positions (the corresponding forces on the Si atoms 
are identically zero by symmetry). The first principles 
results for the six elastic constants of Q;-Pt2Si are given 
The values labeled as "frozen" correspond 



in Table III 



to keeping all of the atoms held fixed at the positions de- 
termined solely from the strain tensor, while the elastic 
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constants labeled "relaxed" were obtained by relaxing the 
strain-induced forces on the Pt atoms. The bulk mod- 
ulus is calculated from the tetragonal elastic constants, 
Bo ~ ^(2cii -I- C33 + 2ci2 -I- 4ci3), and has the same value 
in the frozen and relaxed calculations. 

The requirement that the crystal be stable against any 
homogeneous elastic deformation places restrictions on 
the elastic constants, just as in the cubic case. For tetrag- 
onal crststals these mechanical stability restrictions are as 
foUowsE 



(cii - C12) > 0, (cii + C33 - 2ci3) > 0, 

Cii > 0, C33 > 0, C44 > 0, C66 > 0, 

(2cii -I- C33 -I- 2ci2 + 4ci3) > 0. 



(4) 



These stability conditions again lead to restrictions on 
the magnitude of Bq (Ref. ^ 



1. 

3(^12 



2c 



13, 



<Bo< i(2cii 



C33)- 



(5) 



There are nine independent and non-zero elastic con- 
stants for orthorhombic PtSi. Relaxation of the inter- 
nal degrees of freedom was necessary in calculating all 
nine PtSi elastic constants because the atomic positions 
are not completely fixed by the space group symmetry, 
even for the unstrained crystal, and consequently there 
exist free internal parameters (see Table |) which must 
be redetermined for any distortion of the crystal, includ- 
ing hydrostatic pressure. The results of the calculations 
are listed in Table III. As in the case of a-Pt2Si, the 



values labeled as "frozen" correspond to calculations in 
which the internal structural parameters upt/si a-nd Upt/si 
were held fixed at their self-consistent equilibrium values, 
while in the "relaxed" calculations these internal struc- 
tural parameters were recalculated to minimize the en- 
ergy. The values of Bq are obtained from the elastic 
constants, B^ ^ i(cii -f- C22 + C33 + 2ci2 + 2ci3 + 2C23). 
As expected, the relaxed value of Bq is smaller than the 
frozen value. 

Mechanical stability leads to restrictions on the elas- 
tic coiKjants, which for orthorhombic crystals are as 
followsO 

(cii + C22 - 2ci2) > 0, (cii -t- C33 - 2ci3) > 0, 

(C22 + C33 - 2C23) > 0, 

cii > 0, C22 > 0, C33 > 0, (6) 

C44 > 0, C55 > 0, C66 > 0, 

(cii -I- C22 + C33 -I- 2ci2 -I- 2ci3 + 2C23) > 0. 

As in the case of Q!-Pt2Si, we can obtain restrictions on 
the magnitude of Bq, 



1. 

3(012 



Ci3 



C23) < -Bo < i(cii + C22 + C33). (7) 



From Tables || and |ni we see that overall the elas- 
tic constants of a-Pt2Si and PtSi appear much closer to 
those of pure Pt as compared to pure Si. The trends of 



the elastic constants as a function of the atomic percent 
Pt in all four materials are plotted in Fig. ^. Each of 
the curves corresponds to an average of a different class 
of elastic constants, while the symbols show the values 
of the individual elastic constants themselves. As we see 
from Eqs. (^), (||) and (0), mechanical stability requires 
that Bq be larger than the average of cn, C22, and C33 
but smaller than the average of C12, C13, and C23 [note 
that in the case of Q;-Pt2Si the appropriate averages are 
■i(2cii + C33) and ^{ci2 + 2C13) because Cn = C22 and 
C13 — C23 for tetragonal crystals] . This stability require- 
ment is reflected in the top three curves in Fig. ||. We 
also see that these three curves each increase monoton- 
ically as a function of atomic percent Pt from pure Si 
to pure Pt and we note that all three classes of elas- 
tic constants represented by these curves correspond to 
strains in which the volume is not fixed. Conversely, 
the two lower curves labeled (cn — ci2)/2 and C44 corre- 
spond to the two classes of elastic constants in which the 
strains are strictly volume-conserving [in the case of PtSi 
the lowest solid-line curve and large open circles corre- 
spond to elastic constant combinations -1(011-1-022 — 2012), 
3:(cii-l-C33-2ci3), and 5(022+033-2023)]. We see that in 
this case the two sets of averages are approximately con- 
stant as a function of atomic percent Pt. The significance 
of this difference in the trends of volume-conserving ver- 
sus non- volume-conserving elastic constants is connected 
to the curve labeled Co and is discussed in Sec. ^ along 
with a general discussion of the relationship between the 
magnitudes of the various elastic constants and the chem- 
ical bonding. 



IV. ELECTRONIC CHARGE DENSITY 



In order to provide insight into the nature of the chem- 
ical bonding in a-Pt2Si and PtSi we have analyzed the 
valence electronic charge density in these materials. We 
have chosen to plot charge density differences, the su- 
perposition of free atom densities subtracted from the 
fully self-consistent crystal density, thus emphasizing the 
formation of bonds. Since we are using an all-electron 
method the calculated charge density has a large am- 
plitude close to each of the atomic cores. Subtracting 
two such large numbers can sometimes produce unusual 
features in the plots described below, but these are of 
no consequence to our discussion. Rather we focus on 
the smoothly varying density differences in between the 
atomic cores. In all of the gray-scale plots the brighter 
spots represent an increase in the density relative to su- 
perimposed free atoms while the darker spots represent 
a decrease, with the exact same scale being used in all of 
the plots. 
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A. FPLMTO method 



B. Pt and Si 



The valence electronic charge densities were obtained 
using a fn.ll-potential linear muffin-tin orbital (FPLMTO) 
methodtSilj which makes no shape approximation for the 
crystal potential. The crystal is divided up into regions 
inside atomic spheres, where Schrodinger's equation is 
solved numerically, and an interstitial region. The wave- 
functions in the interstitial are Hankel functions. An 
interpolation procedure is used for evaluating interstitial 
integrals involving products of Hankel functions. The 
triple-K basis is composed of three sets of s, p, d, and 
/ LMTOs per atom with Hankel function kinetic ener- 
gies of -K^ = -0.01, -1.0, and -2.3 Ry (48 orbitals 
per atom). The Hankel functions decay exponentially 
as e~'^^. The angular momentum sums involved in the 
interpolation procedure are carried up to a maximum 
of ^ = 6. The calculations presented here are based 
on the local density approximation (LDA), using the 
exchange-correlation potential of Ceperley and_AlderE2l 
as parameterized by Vosko, Wilk, and Nusair.tJ The 
scalar-relativistic Schrodinger equation was solved self- 
consistently. We did not include spin-orbit interactions 
and we used atomic sphere radii equal to one-half the 
nearest-neighbor bond lengths. In the case of a-Pt2Si 
we included an empty atomic sphere at the octahedral 
interstitial site, as well as the usual empty spheres in the 
interstitial sites of pure cubic-diamond-phase Si. How- 
ever, these empty spheres do not contribute to the basis 
but merely improve the accuracy of the interstitial inter- 
polation procedure. 

The Pt 6s, 6p, 5d, and 5/ orbitals as well as the Si 
3s, 3p, 3d, and 4/ orbitals were all treated as valence 
states. The semi-core Pt 5s and 5p orbitals were treated 
as full band states by carrying out a "two-panel" cal- 
culation. The second panel band calculation for the 
semi-core orbitals included the Pt 5s, 5p, 5d, and 5/ 
orbitals as well as all of the Si valence orbitals. The Bril- 
louin zone (BZI sums were carried out using the tetrahe- 
dron method.tJ We used the same mesh of k-points for 
both the self-consistent total energy and charge density 
calculations. In the case of Q!-Pt2Si we used a shifted 
24x24x24 (12x12x12) mesh in the full BZ, resulting in 
1056 (159) irreducible k-points in the first (second) panel. 
In the PtSi calculations we used a shifted 12x16x12 
(6x8x6) mesh in the fuU BZ, resulting in 288 (36) ir- 
reducible k-points in the first (second) panel. A shifted 
28x28x28 (16x16x16) mesh in the full BZ was used for 
fee Pt, resulting in 2030 (408) irreducible k-points in the 
first (second) panel. Finally, a shifted 12x12x12 mesh 
in the full BZ was used for cubic-diamond-phase Si, re- 
sulting in 182 irreducible k-points. 



We start with the well-known cases of pure diamond- 
phase Si and fee Pt in order to provide a baseline with 
which to compare the results we obtain for the silicides. 
In Fig. ^(a) we see the localized piling up of additional 
charge between each pair of Si atoms which corresponds 
to the covalent bonds in this material. Except for these 
bonds, the density is relatively unchanged from the free- 
atom superposition in the remaining regions outside of 
the atomic cores, as can be seen by identifying the "0" 
level in the accompanying scale bar. This circumstance 
is in stark contrast to the case of fee Pt in Fig. ^(b). In 
Pt the increase in density is spread approximately uni- 
formly throughout all of the regions outside the atomic 
cores. In fact, from this perspective Pt appears almost 
free-electron-like, despite the more localized nature of the 
states arising from the partially occupied d-band. Thus 
we see that charge density difference plots such as those 
in Fig. ^ are clearly able to distinguish metallic bond- 
ing, as occurs in Pt, from covalent bonding, as occurs in 
Si. For later purposes we note that the nearest-neighbor 
spacing is 2.35 A in Si and 2.77 A in Pt. 



C. a-PtaSi 

Each Si atom in a-Pt2Si has eight Pt nearest-neighbors 
at a distance of 2.47 A [see Fig. |l|(a)]. In addi- 
tion to four Si nearest-neighbors, each Pt atom also 
has four Pt second-nearest-neighbors at a distance of 
2.79 A and two Pt third-nearest-neighbors at 2.98 A. The 
Pt second-nearest-neighbors form two-dimensional (001) 
planes while the Pt third-nearest-neighbors form linear 
[001] chains. The Pt second-nearest-neighbor distance is 
very close to the nearest-neighbor distance in pure fee Pt 
and thus we might expect these two-dimensional planes 
to exhibit evidence of metallic bonding. This is in fact 
what we see, as shown in Fig. ^(a) which bears a strong 
resemblance to the analogous plot in Fig. §(b). However, 
this approximately uniform increase in the charge den- 
sity in the regions outside the atomic cores is confined 
to the two-dimensional second-nearest-neighbor Pt (001) 
planes. In particular, there is little evidence of bonding 
(i.e. little or no increase in the charge density relative 
to free atoms) along the third-nearest-neighbor [001] Pt 
chains. 

In addition to the two-dimensional "metallic" bond- 
ing, we find strong evidence of covalent bonding between 
the Pt and Si nearest-neighbors, illustrated in Fig. ^(b). 
Unlike the case of pure Si where the increase in charge 
density occurred between pairs of atoms, here the density 
increase is localized between three atoms, two Pt and a 
Si. For this reason we refer to these as three-center co- 
valent bonds. We might even be tempted to call these 
four-center bonds because there is a smaller increase in 
the density, in between the two Pt atoms, which connects 
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two of the three-center bonds. However, we note that 
the X-axis in Fig. ^(b) is along the [110] direction and 
that each of the Pt-Pt pairs in between two of the three- 
center bonds are also located in one of the (001) planes 
which exhibit evidence of metallic bonding [Fig. H(a)]. It 
would thus appear that rather than four-center covalent 
bonds, a more appropriate description of the bonding in 
a-Pt2Si would be three-center bonds interconnected by 
two-dimensional metallic sheets. 

The y-axis of Fig. ^(b) is along the [111] direction 
which highlights two of the central Si atom's three-center 
bonds. However, from Fig. |l|(a) we see that there are 
four of these crystallographic directions and therefore a 
total of 8 of these three-center bonds for each Si atom. 
As noted above, the pair of Pt atoms participating in 
a given three-center bond are second-nearest-neighbors 
themselves. Fig. ^(c) shows that there is another set of 
three-center bonds involving one Si atom and a pair of Pt 
atoms which are third-nearest-neighbors oriented along 
the [001] chains. The x-axis in Fig. |(c) is along [100] 
and the y-axis is along [001]. There is little or no indi- 
cation of an increase in charge density along the Pt-Pt 
[001] chains. In addition to the two three-center bonds 
in Fig. ^(c), there are two more of these bonds located in 
the plane obtained by a 90° rotation about the [001] axis 
see Fig. 0(a)], for a total of four of these three-center 
bonds for each Si atom. 

Thus wc sec that each Si atom in a-Pt2Si participates 
in twelve three-center bonds, eight with Pt-Pt second- 
nearest-neighbors and four with Pt-Pt third-nearest- 
neighbors, and that these three-center covalent bonds 
are interconnected by two-dimensional second-nearest- 
neighbor Pt-Pt metallic sheets. Given the large increase 
in the number of bonds in a-Pt2Si relative to pure Si 
we expect that each individual bond will be weaker than 
one of the covalent bonds in Si. However, taken as a 
whole and in terms of the material strength, the more 
distributed nature of the bonding in a-Pt2Si may indi- 
cate something closer in character to the pure metallic 
bonding in fee Pt. This interpretation is supported by 
the calculated elastic constants in Fig. ||, where the non- 
volume-conserving elastic constants for a-Pt2Si are much 
closer to those of fee Pt as opposed to pure Si. We ad- 
dress this issue in more detail in Sec. 0. 



D. PtSi 



In the orthorhombic PtSi structure each Si atom has 
six Pt neighbors, with one Pt at 2.41 A, two at 2.43 A, 
one at 2.52 A, and two at 2.64 A. In view of the fact 
that the nearest-neighbor Pt-Si distance is 2.47 A in a- 
Pt2Si it is perhaps not surprising that we find the two 
Pt neighbors at 2.64 A appear to contribute little to the 
bonding in PtSi. Each Si also has two Si fifth-nearest- 
neighbors at 2.84 A but again we find little evidence of 
bonding between these atoms which is consistent with 



the fact that the nearest-neighbor distance in pure Si 
is only 2.35 A. In addition to six Si neighbors at the 
same distances listed above, each Pt atom also has two Pt 
neighbors at a sixth-nearest-neighbor distance of 2.87 A 
and two more at a seventh-nearest-neighbor distance of 
2.90 A. These distances are somewhat larger than the 
2.77 A nearest-neighbor distance in pure fee Pt. 

The striking appearance of three-center bonds in a- 
Pt2Si is repeated in orthorhombic PtSi, as shown in 
Fig. ||(a). As we see in Fig. |l](b), a convenient way to 
think of the PtSi structure is as two alternating planes 
of atoms stacked along the b-axis. Fig. ^(a) shows the 
charge density difference in one of these planes. As in the 
case of a-Pt2Si [Figs, ^(b) and ^(c)] we see a pileup of 
charge relative to the free atom density which is not local- 
ized between a single pair of atoms but rather between 
one Si and two Pt atoms. These Pt neighbors partici- 
pating in the three-center bond are the first- and third- 
nearest-neighbors of the Si atom and are at distances of 
2.41 A and 2.52 A. The two Pt atoms are themselves 
sixth-nearest- neighbors, with a bond length of 2.87 A. 
There appears to be a small increase in the charge den- 
sity between these two Pt atoms. We note that the two 
different three-center bonds shown in Fig. |^(a) are equiv- 
alent by symmetry. 

The two second-neighbor Pt atoms of a given Si atom 
are located in adjacent b-axis planes from the Si. The 
charge density difference for these bonds is shown in 
Fig. |5|(b), which indicates that they are of the stan- 
dard two-center variety. In addition to these two-center 
bonds, the plot also shows part of the bond with the first- 
neighbor Pt atom on the left side of the figure. In fact, 
the two second-neighbors as well as the first- and third- 
neighbors form a very distorted tetrahedron around the 
central Si atom. The Pt-Si-Pt bond angles involving 
one Pt second-neighbor and one third-neighbor are very 
nearly equal to the perfect tetrahedral angle of 109.47° 
in pure Si, but the remaining four bond angles vary con- 
siderably, ranging from 71° to 131°. 

There is very little evidence of an appreciable increase 
in the charge density between the Pt and Si fourth- 
nearest- neighbors and the Si-Si fifth-neighbors, as we 
mentioned above. The Pt-Pt sixth-nearest-neighbors in 
Fig. ||(a) show some evidence of charge accumulation 
but the Pt-Pt seventh-neighbors do not. We thus see 
that there appears to be only two sets of strong covalent 
bonds in orthorhombic PtSi, the three-center Pt-Si-Pt 
bonds within a given b-axis plane and the two-center 
Pt-Si bonds between atoms in adjacent b-axis planes, 
resulting in a total of only three bonds per Si atom. In 
this sense the bonding in PtSi appears to be qualitatively 
much more similar to that in pure Si as compared to pure 
Pt or even Q;-Pt2Si. In particular, we are unable to iden- 
tify any concrete evidence in PtSi of a uniform increase in 
interstitial charge density that might be associated with 
an element of metallic bonding. We revisit this subject 



in Sec. VD 
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V. VALENCE FORCE FIELD MODELS 



C44, 



In order to provide a more quantitative analysis of the 
trends in the elastic constants as well as the various ele- 
ments of the chemical beading , we construct simple va- 
lence force field modelscaa to describe the interatomic 
interactions for pure Pt, pure Si, and the two silicides. 
In these models the change in energy upon distorting the 
crystal, Ai?, is given as follows 



1 



N 



2 

N 



(8) 



where N is the number of primitive cells in the crystal, n 
is the number of atoms in the primitive cell, V is the vol- 
ume, Ac?i is the change in the ith bond length , and A0y 
is the change in the bond angle between the zth and jth 
bonds. We determine the C and K force constants by 
equating this expression for Ai? to corresponding elas- 
tic constant expressions derived from Eq. (|l|), examples 
of which are given in Ref. |^. These coefficients are re- 
ferred to as force constants because Eq. (^) could also be 
used to anaylze the phonon spectrum and in this case, 
within a constant factor, the coefficients play the role of 
Hooke's law force constants. The factor of n in the first 
term of Eq. (||) is explicitly included so that the result- 
ing force constant Co represents the volume contribution 
per atom, thus facilitating the comparison between ma- 
terials with different numbers of atoms in the primitive 
cell. Similarly, the indices i and j are summed over all 
of the relevant bonds for each of the atoms in the prim- 
itive cell (avoiding any double counting), which results 
in force constants that represent the interaction strength 
for a single bond (C^) or bond angle (Kij). The vol- 
ume term in Eq. (^ is needed for metals such-as-Pt and 
is reminiscent of the embedded-atom methodOES which 
has been successful in treating fee metals. Similarly, the 
angular terms are needed for covalently bonded systems 
such as Si; sucb-J.erms are a part of the Tersoff poten- 
tial formulationEZi which has been used successfully in 
semiconductor systems. Both the volume and the angiu 
lar terms lead to deviations from the Cauchy relationsEa 
which are strict equalities between various elastic con- 
stants that apply when the interatomic interactions are 
purely pairwise [i.e. including only the second term in 
Eq. (D]. 



A. Pt 

In the case of fee Pt we construct a two-parameter 
model, considering only the nearest-neighbor bond length 
and a volume term, but no angular terms. The ra- 
dial force constant Ci can be obtained from the volume- 
conserving strains corresponding to either (cn — C12) or 



(cii - C12) 



1 1 

V 4 



Ci, 



and 



where 



C44 



1 1 

V 2 



Ci, 



(9) 



(10) 



1„3 



a is the volume per atom. Taken together, 
these two equations provide an explanation for the fact 
that (cii — C12) and C44 for pure Pt are similar in mag- 
nitude in Table ||. They also satisfy the cubic stability 
requirements [Eq. (^)] that (cn — C12) > and C44 > 0. 

In order to facilitate comparison with the silicides 
where there is no experimental data for the elastic con- 
stants, we also use the theoretically determined elastic 
constants for fee Pt. For the purpose of internal consis- 
tency we use the theoretical equilibrium volume as well. 
The resulting two values of Ci obtained from Eqs. (|^) 
and ( p^ ) are 15.78 eV and 16.42 eV, respectively. The 
fact that the two numbers differ is an indication of the 
incompleteness of the two-parameter model. Use of the 
experimental elastic constants yields a bigger difference 
but we nonetheless will use the average of these two val- 
ues for the purpose of comparing to the silicides. 



Ci =v [(cii - C12) -I- C44] 



(11) 



Evaluated using the theoretical elastic constants we ob- 
tain Ci — 16.10 eV, while the experimental elastic con- 
stants correspond to a value of 16.95 eV. We note that 
we could have eliminated the need to use the averaged 
expression in Eq. ( pT| ) by including additional force con- 
stants but we prefer to maintain the conceptual simplicity 
of the two-parameter model. For example, including an 
angular interaction in Eqs. ^ and (|l^ results in a small 
and slightly negative angular force constant K which is 
conceptually unsatisfying. 

The uniform expansion and compression represented 
by the bulk modulus Bq can be used to obtain the fol- 
lowing expression involving Cq and Ci , 



2ci 



Co + ttCi 



(12) 



Eqs. (|ri|) and (^ together yield a value of Co = 16.54 eV 
using the theoretical elastic constants, and 15.73 eV us- 
ing the experimental values. The values of Co and Ci 
obtained from the theoretical elastic constants are listed 
in Table IV . We note that the volume force constant has 



approximately the same magnitude as the radial force 
constant and that both are important in contributing to 
the large bulk modulus. For completeness we also give 
the expressions for cn and C12 in terms of Co and Ci, 



cii = -(Co-I-Ci), 

V 



(13) 
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Cl2 



Co + Ici 



(14) 



Eqs. (|12|-|14D explicitly satisfy the stability requirement 
that ci2 <Bo< cii [Eq. (1)]. 

Tiie Cauchy relation for cubic crystals is that ci2 — 
C44B3 Using Eqs. ( [To| ) and (|l^ we obtain the follow- 
ing expression for the deviation from the cubic Cauchy 
relation, 



The final independent elastic constant is the bulk mod- 
ulus Bq which corresponds to an isotropic expansion or 
compression and therefore only involves radial but not 
angular distortions. In addition, this distortion is not 
volume-conserving and thus wc could also have included 
the Co volume term from Eq. (^) , which would not affect 
either of the volume-conserving strains corresponding to 
Eqs. (16) and (pj|), but would yield the following equa- 
tion for the Bq distortion, 



(ci2 - C44) = -Cq. 

V 



(15) 



Thus we see that the large and positive deviation from 
the Cauchy relation in pure Pt is due to a large volume 
contribution to C12. Moreover, the presence of the vol- 
ume contribution Cq is responsible for the fact that Cn, 
C12, and Bq are all significantly larger than the volume- 
conserving elastic constants ^(cn — C12) and C44. 



B. Si 

In the case of Si we also construct a two-parameter 
model but instead consider only the nearest-neighbor 
bond length and the tetrahedral bond angle and set all 
of the other C and K force constants to zero. Since 
there are three elastic constants and we allow only two 
force constants, we can check the accuracy of the model. 
The volume-conserving strain corresponding to (cn — C12) 
leaves the nearest-neighbor bond lengths unchanged to 
first order in the distortion and thus only the angular 
force constant enters. 



1 



(cii 



C12) = - 2Kn, 

V 



(16) 



where w = is the volume per atom. Comparing this 
result to Eq. (||) we see that ^(cn — C12) has a very dif- 
ferent origin in Si as compared to fee Pt, despite the fact 
that the two values are approximately the same in Ta- 
ble g 

The volume-conserving C44 strain in Si involves both 
radial and angular distortions and can thus be used in 
conjunction with Eq. (Hq) to determine Ci, 



frozen 
C44 



-Ci + -Kn 



(17) 



where c^™^''" corresponds to a pure C44 strain, without 
allowing for any internal relaxation. This choice is con- 
venient but not essential and we can test how well the 
two force constants describe the final remaining elastic 
constant. Using the theoretically determined elastic con- 
stants and equilibrium volume, Eqs. (|l^ and ( p7| ) yield 
Ci = 54.06 eV and Ku = 3.13 eV, which are listed in 
Table Based on this analysis, the fact thatfcn — C12) 
and C44 for Si are close in magnitude in Table O is merely 
a coincidence having to do with the specific values of the 
Ci and Ku force constants. 



Bo = l{cn+2c,2) = -(co + lc, 
3 V \ 9 



(18) 



Our two-parameter model has Cq = and thus the extent 
to which Co obtained from Eq. ( p^ ) deviates from zero 
provides a direct measure of how well the two-parameter 
model is able to describe the elastic constants. Using the 
theoretical values determined here, Eq. (|l^) yields Co = 
—0.18 eV, which demonstrates that the two-parameter 
model is indeed sufficiently accurate for describing the 
elastic constants in Si. For the sake of completeness we 
give the expressions for cn, C12 and the deviation from 
the Cauchy relation (ci2 — C44), including a volume con- 



tribution, 



Cll 



C12 



Co + -Ci + -A'l] 



Co 



2 4 
9^^ - 3^" 



(Ci2 - C44) 



frozen f ^^7-^ 

V \ 9 



(19) 



(20) 



(21) 



We can compare 0|Hif two-parameter model to the one 
derived by Harrison.E^ His angular term has the iden- 
tical form as ours and his value of Ku = 3.2 eV differs 
from ours of 3.13 eV only because we have used the theo- 
retical elastic constants and equilibriuui lattice constant 
while he uses the experimental values.LJ We derived the 
value of the radial force constant Ci using the frozen 
C44 elastic constant whereas Harrison derives his radial 
force constant from Bq. The two values would be identi- 
cal if the value of Co derived from Eq. ( p^ were exactly 
zero. The small deviation from zero, in addition to the 
difference in the lattice constants used, leads to a small 
difference between Harrison's value of Ci — 55.0 eV and 
our value of 54.06 eV. 

We have found that for Si the angular force constant 
Ku is more than an order of magnitude smaller than 
the radial force constant Ci (Table IV ) . The angular in- 
teraction is nonetheless of particular importance for two 
reasons. The first is that the crystal would be unstable in 
the absence of angular interactions since (cn — C12) = 
[Eq. 1^)] and Bq = cn = C12 [Eqs. <^^^] for purely 
radial interactions, both of which violate the stability 
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conditions discussed in Sec. III. The second reason is 



that the angular interaction is responsible for the fact 
that the elastic constants do not obey the Cauchy rela- 
tion for a cubic crystal, C12 — C44 [Eq. Including 
a volume term but not an angular interaction would still 
result in an unstable crystal since the (cn — C12) distor- 
tion is volume-conserving and the dependence on Co is 
identical for Bq, cn, and C12. In addition, a positive 
volume contribution in the absence of an angular inter- 
action could not account for the fact that the deviation 
from the Cauchy relation is negative (ci2 < C44). Thus 
we see that Eq. ( ^l| ) and the fact that Ku is small com- 
bine to provide an explanation for the fact that C12 is 
smaller than C44, but only by a relatively small amount. 
This is in sharp contrast to pure Pt where the deviation 
from the Cauchy relation [Eq. ([Tsl)] arises from Co and 
is large and positive, resulting in a value of C12 which is 
more than four times larger than in Si. More generally, 
the absence of a volume contribution in Si is responsible 
for the fact that C12 and Bq are similar in magnitude to 
the volume-conserving elastic constants ^(cii — C12) and 
C44, in contrast to the case of Pt. 

From the force constants listed in Table |^ we see that 
Ci for Si is more than a factor of three times larger than 
for fee Pt which is consistent with the presence of strong 
covalent bonds in Si and distributed metallic bonding in 
Pt. In addition, the volume contribution Co is equally 
important in terms of the metallic bonding in Pt but 
plays no role in Si. The influence on the elastic constants 
of these qualitative differences in the chemical bonding 
are clearly illustrated by comparing Eqs. (|l^) and ( p^ ) 
for the bulk moduli in Pt and Si, respectively. We see 
that the geometry coeflScient of Ci is three times larger 
for Pt than for Si, reflecting the difference in the nearest- 
neighbor coordination and nearly compensating for the 
difference in the magnitudes of the two force constants. 
Given that Co and Ci are approximately the same in Pt, 
we see from Eq. (jl^) that the volume contribution to Bq 
is approximately 50% larger than the contribution from 
Ci. In the case of Si the Co force constant is essentially 
zero and this difference accounts for most of the differ- 
ence in the magnitudes of Bq between Pt and Si. The 
pre- factors of 1/v account for the remaining difference 
since the volume per atom v is 30% larger in Si. 

We therefore see that the presence or absence of metal- 
lic bonding, as reflected in the Co and Ci force constants, 
is intimately connected to the magnitudes of Bq. Simi- 
lar analyses can be used to explain the fact that cn and 
C12 are also larger in Pt, the predominant reason being 
the presence of a large volume contribution (or equiva- 
lently, metallic bonding) in Pt but not Si. Conversely, the 
elastic constants corresponding to both of the volume- 
conserving distortions in Pt, ^(cn — C12) and C44, are 
approximately the same as in Si, indicating that they 
are less sensitive to the differences in chemical bonding 
for these two materials. These differing trends in the 
volume-conserving versus non-volume-conserving elastic 
constants were already noted in Sec. [II and are illus- 



trated in Fig. ||. The volume force constant Co is in- 
cluded in the figure on the same scale as the elastic con- 
stants by dividing by the appropriate volume per atom 
V (note that -Cq is precisely the combination that en- 
ters all of the expressions for the non- volume-conserving 
elastic constants). 



C. a-PtaSi 



In order to provide a more quantitative description of 
the chemical bonding in a-Pt2Si we describe the inter- 
atomic interactions using a valence force field model, just 
as we did for pure fee Pt and pure Si. In view of the 



analysis of the valence charge density in Sec. IV C , we 
include first-, second-, and third-nearest-neighbor radial 
force constants as well as a volume term. We also con- 
sider some of the angular interactions. In keeping with 
our neglect of angular interactions in pure Pt we also 
neglect the bond angles between any two Pt-Pt bonds, 
both in the two-dimensional (001) metallic sheets and 
the [001] Pt chains. In consideration of the three-center 
bonds discussed earlier, we include both the Pt-Si-Pt 
and the Si-Pt-Pt bond angles relevant to the three-center 
bonds involving one Si atom and two second-neighbor 
Pt atoms. However, we neglect the bond angles rele- 
vant to the three-center bonds involving one Si and two 
third-neighbor Pt atoms. This choice is based on the ex- 
pectation that the strength of the angular interactions 
will generally be smaller than that of the radial inter- 
actions and that the three-center bond involving two Pt 
second-neighbors is stronger than the one involving two 
third- neighbors. We thus have six force constants which 
can be fit to the six elastic constants. For the sake of 
convenience we fit the force constant expressions to the 
"frozen" elastic constants, where no internal relaxations 
were carried out. This choice is not essential and need 
not be considered an additional approximation because 
the resulting force constants could be used to directly 
calculate the internal relaxations. 

The volume-conserving strains corresponding to (cn — 
C12) and C44 both depend only on the first-nearest- 
neighbor Pt-Si radial force constant Ci since the second- 
and third-nearest-neighbor Pt-Pt bond lengths are left 
unchanged to first order. In addition, (cn — C12) depends 
on the Si~Pt-Pt bond angle but not the Pt-Si-Pt bond 
angle, while C44 depends on both. We label the force 
constant for the Pt-Si-Pt bond angle as Ku because it 
is the angle between two first-neighbor bonds. Similarly 
we label the Si-Pt-Pt force constant as i^i2. Equating 
the elastic constant [Eq. (^] and force constant [Eq. (H)] 
expressions for the change in the energy and using the 
theoretical lattice constants from Table ||, we obtain the 
following two equations, 

lien - ci2)f"^^" = -(0.2685 Ci + 1.1553 K,^), (22) 
2 V 
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and 



frozen ^ _(o.3092 Ci + 0.2874 Ku + 0.8089 K12), 

V 



(23) 



where v = 0.2529 is the volume per atom. We note 
that in these equations and aU of those that foUow, the 
numerical coefficients are simply geometrical factors con- 
taining various combinations of the a and c lattice con- 
stants. Examining the geometry coefficients of Ci in 
Eqs. and (p^), together with the expectation that 
the angular force constants will be significantly smaller 
in magnitude than Ci, we see that these two equations 
provide a natural explanation for why ^(cn — C12) is 
si milar in magnitude but slightly smaller than C44 in Ta- 
ble We also see that two of the stability requirements 
[Eq. (y)], (cii — C12) > and C44 > 0, are explicitly sat- 
isfied. 

The volume-conserving strain corresponding to cge 
changes the second-neighbor Pt-Pt bond length, leav- 
ing the other two bond lengths unchanged to first order. 
This strain also modifies the two bond angles yielding, 

= i + 1.2394 Kii + 0.6197 Ki^j ■ (24) 



Thus we can see from Eqs. (|2^) and (^^ that since cge in 
Table [I] is somewhat more than two times larger than 
C44, we expect that the second- neighbor Pt-Pt force con- 
stant C2 must be approximately two times larger than 
the first-neighbor Pt-Si force constant Ci. We will in 
fact find this to be the case. Eq. ( |24| ) also satisfies the 
stability requirement that cqq > 0. 

The final remaining volume-conserving strain corre- 
sponding to (cii + C33 — 2ci3) changes all of the first-, 
second-, and third-neighbor bond lengths, as well as the 
two bond angles, 



^(cii + C33 - 2ci3) 



frozen 



i f 0.0687 Ci + -^C2 + ic3 
V \ 12 6 



+ 0.3723 Kn + 0.4750 K 



12 



(25) 



We note that j{cii + C33 — 2ci3) > as required for me- 
chanical stability, and that it is similar in magnitude to 
i(cii — C12) and C44 in Table 

The uniform expansion and compression correspond- 
ing to the bulk modulus Bq changes the volume and all 
of the bond lengths but leaves the bond angles fixed. 



Bo = -(2cii + C33 + 2ci2 + 4ci3) 



(26) 



The final two equations resulting from the cn and C33 
strains both include a contribution from a change in the 
volume, 



1 



rrozcn ^ _ I (J^ ^ 0.5370 Cl 

V 



1 



and 



-C2 + 0.1655 Kn + 1.2381 



C33 = -(co + 0.3560 Cl + 
V \ 3 



(27) 



0.6619 Xii +0.3310 Ki2 



Similarly, the equations for ci2 and C13 are 



(28) 



1 / „ 1 
.V 3^ 



^frozen ^i^(jo + -C2 + 0.1655 Kn - 1.0726 K12 ] , 



(29) 



and 



1 



ci3 = -(Co + 0.3092 Cl - 0.3310 ^"11 - 0.1655 fCia). 

V 

(30) 

Eqs. (|26|)-(|3^ explicitly satisfy the stability require- 
ments [Eq.^)] that Bq < ^(2cii -|- C33) and Bq > 
i(ci2 -I- 2C13). We note that in contrast to the case of 
Si, all of the stability requirements in Eq. (^ would be 
satisfied even for purely radial interactions (i.e. no angu- 
lar interactions). 

Eqs. (|22|)-(|2|), (|^, and (||) represent six hnearly in- 
dependent equations in the six unknown force constants. 
Solving this linear system of equations yields the force 
constants listed in Table |l^ for Q;-Pt2Si. The volume 
force constant Cq is only 16% smaller than in pure Pt. 
This is consistent with the presence of two-dimensional 
metallic sheets in a-Pt2Si and the fact that there are 
a large number of distributed three-center bonds all in- 
terconnected by these sheets. The first neighbor Pt-Si 
force constant Ci is nearly four times smaller than Ci in 
pure Si. This large reduction results from the fact that 
each Si atom in a-Pt2Si has eight Pt nearest-neighbors 
and participates in twelve different three-center bonds. 
Conversely, the second-neighbor Pt-Pt force constant C2 
is 60% larger than the corresponding Ci force constant 
in pure fee Pt, despite the fact that the two Pt-Pt bond 
lengths are very nearly the same. We can understand this 
result because each Pt atom in pure Pt has twelve near- 
est neighbors while each Pt in Q!-Pt2Si has only four Pt 
second-neighbors and two Pt third-neighbors. Moreover, 
the Pt atoms in the silicide participate in covalent three- 
center bonds in addition to the metallic bonding within 
the two-dimensional sheets. The distributed nature of 
these bonds and the large number of them in the primi- 
tive cell are both consistent with the fact that C2 is still a 
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factor of two smaller than Ci in pure Si. We found little 
evidence of an increase in the electronic charge density 
between the Pt-Pt third-neighbors and this is reflected 
in the fact that C3 is more than three times smaller than 
C2. We also find that the angular force constants Ku 
and K12 are similar in magnitude to the Ku force con- 
stant in pure Si. These angular terms play an important 
but less crucial role in the silicide as compared to pure 
Si. 

Having determined the values of the individual force 
constants we can now use them to understand the trends 
in the elastic constants. For example, the two Cauchy 
relations for tetragonal crystals are that C12 = cge and 
C13 = C44.Ea Using Eqs. (^J) and ( p9| ) the deviation from 
the first Cauchy relation is given by 



(ci2 - cee) 



frozen 



1 



(Co - 1.0739 Kn - 1.6923 K^) 



(31) 



Similarly, the deviation from the second Cauchy relation 
is 



(Cl3 - C44) 



frozen 



1 



(Co - 0.6184 Kn - 0.9744 K^). 



(32) 



As in the case of pure Pt [Eq. (p^)], it is the presence of 
the volume interaction which produces a positive devia- 
tion from the Cauchy relations. The angular interactions 
provide a negative contribution to Eqs. (|T]) and (32), 
just as they did for pure Si [Eq. (pT|)]. From the geome- 
try coefficients of the angular terms we see that (013 — C44) 
must be larger in magnitude than (ci2 — cge). The net re- 
sult is that the deviations from the Cauchy relations for 
a-Pt2Si are still positive but are factors of 2-3 smaller 
than the deviation in pure Pt. 

In pure Pt we found that the volume- conserving elas- 
tic constants were all significantly smaller than the others 
and that this was due predominantly to the presence of 
a large volume contribution Cq. We find the same trend 
in a-Pt2Si with i(cii - C12), C44, and i(cii -I- C33 - 2ci3) 
all being similar in magnitude and smaller than all the 
remaining elastic constants (see Fig. |]). The notable ex- 
ception to this trend is cee- In conjunction with Eq. ( |2^ ) 
we already noted that the large value of cee in relation to 
the other volume-conserving elastic constants is due pri- 
marily to the fact that the second-neighbor Pt-Pt force 
constant C2 is a factor of two larger than the first neigh- 
bor Pt-Si force constant Ci. This result is in turn di- 
rectly related to the presence of the network of three- 
center bonds interconnected by two-dimensional metallic 
sheets. We also saw that mechanical stability requires 



C12 



< Cii, Ci3 < ^(cii 



C33) 



Bo < i(2cii + C33), and 



Bo > ^{ci2 + 2C13) and that the deviations from the 
Cauchy relations are positive. The remaining variations 
among the six elastic constants in Table [II are deter- 
mined by the detailed dependence on the various force 
constants as described above. 



One interesting example is that C13 is essentially iden- 
tical to C66. Comparing Eqs. ( ^ and (|30| ) we see that in 
the case of C13 positive volume and first-neighbor radial 
terms are partially counterbalanced by negative angu- 
lar contributions whereas cge corresponds to a volume- 
conserving distortion and has positive angular contribu- 
tions. In addition, the C13 distortion changes the first- 
neighbor bond lengths, leaving the others fixed, while the 
C66 distortion changes the second- neighbor bond lengths, 
leaving the others fixed. We have already noted that cee 
is anomalously large in comparison to C44 predominantly 
because C2 is twice as large as Ci. The elastic constant 
C12 is larger than C13 for the same reason, thus explain- 
ing how it is at least possible for the volume-conserving 
elastic constant cge to be similar in magnitude to the 
non- volume-conserving elastic constant C13 . In summary, 
while we are able to explain the overall magnitudes of 
the individual elastic constants, we are forced to con- 
clude that the specific equality of C13 and cge in Table III 
depends on the precise values of the individual force con- 
stants and is therefore simply accidental. 

We saw in Eq. (|l|) for the bulk modulus of Pt that the 
volume contribution represented 60% of the total, with 
the contribution from the radial interaction making up 
the rest. The same approximate 60:40 split between the 
volume and radial contributions applies to the expression 
for the bulk modulus of a-Pt2Si in Eq. (|2^). In addition, 
the volume per atom v is nearly the same in the two ma- 
terials. Thus we see that the 16% reduction in Co for 
a-Pt2Si relative to Pt, combined with a similar reduc- 
tion in the overall radial contribution, leads to a bulk 
modulus which is approximately 20% smaller in a-Pt2Si. 
As we noted previously, there is no volume contribution 
in Si where the bulk modulus is a factor of 2-3 smaller. 
Conversely, the volume-conserving elastic constant C44 is 
similar in mag nitude for all three materials. Comparing 
Eqs. (10), (|l7|), and (|23| ) we see that in pure Pt C44 arises 



solely from Ci whereas in pure Si and Q;-Pt2Si it arises 
from a combination of Ci and angular contributions. In 
Si the split is 90:10, radial to angular, while in the sili- 
cide the split is only 60:40 since Ci is a factor of four 
smaller. In addition, the volume per atom is 30% larger 
in Si than in the other two materials. The remaining 
elastic constants for a-Pt2Si can be similarly analyzed in 
relation to those of pure Pt and pure Si. 



D. PtSi 



Once again we construct a valence force field model to 
describe the chemical bonding and elastic constant trends 
in PtSi. In keeping w ith th e discussion of the valence 
charge density in Sec. [VD, we include first-, second-, 
and third-neighbor Pt-Si radial force constants (labeled 
Ci, C2, and C3) as well as sixth- and seventh- neighbor 
Pt-Pt radial force constants (labeled Cq and C7). We 
also include a volume term (Co) and three Pt-Si-Pt an- 
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gular force constants. The angular force constants are 
labeled i^ia, corresponding to the bond angle between 
first- and third-neighbor Pt-Si bonds, K22, correspond- 
ing to the bond angle between two second-neighbor Pt-Si 
bonds, and K2'i, corresponding to the bond angle be- 
tween second- and third- neighbor Pt-Si bonds. These 
three bond angles are the ones we have found to be 
most important and are the ones which correspond to 
the distor ted tetrahedral Pt-Si-Pt angles described in 
Sec. [VP. The fourth and last of these angles is rep- 



resented by the force constant K12 but we found it to 
be unimportant and have not included it in the analy- 
sis presented here. Part of the reason for this finding 
may be that this bond angle is 131.72° which is quite 
different from the perfect tetrahedral angle of 109.47°. 
We thus have nine force constants which can be fit to 
the nine elastic constants. As in the case of Q;-Pt2Si we 
fit the force constant expressions to the "frozen" elastic 
constants out of convenience, but this choice is not es- 
sential because the relaxations could be calculated from 
the resulting model. 

Most of the expressions for the elastic constants in 
terms of the force constants involve all of the radial and 
angular terms and thus there is not much to be learned 
by writing them down. Two exceptions are the volume- 
conserving strains corresponding to C44 and cgg which 
depend only on the second-neighbor Pt-Si and seventh- 
neighbor Pt-Pt radial force constants, as well as the an- 
gular force constant i^23- Using the theoretically deter- 
mined structural parameters from Table | we obtain the 
following two expressions, 



frozen 

L44 



1 



-(0.1601 C2 + 0.1182 C7 -f 0.3279 K22,), (33) 



and 



frozen 
^66 



-(0.0882 C2 + 0.00002 C7 -I- 0.1750 



23) 



(34) 



where v = 0.08464 is the volume per atom. The force 
constant C7 will turn out to be small and thus we can 
see from Eqs. (|3^ ) and (HJ) that 044^°'^ is approximately 
a factor of two larger than c|™^°" purely because of ge- 
ometrical factors. In addition, the mechanical stability 
requirements that C44 > and cgg > are satisfied by 
a combination of radial and angular terms. However, 
as in the case of Q!-Pt2Si, the angular terms are not es- 
sential with regard to stability since the crystal would 
still be stable under purely radial interactions. It turns 
out that this circumstance is true for all of the stabil- 
ity requirements in Eq. (H) . We also note that all of the 
volume-conserving elastic constants, ^(cn -I- C22 — 2ci2), 
3:(cii -I- C33 - 2ci3), lic22 + C33 - 2C23), C44, C55, and Cge 
are similar in magnitude and smaller than the other non- 
volume-conserving elastic constants (see Fig. ||). The 
primary exception is c^™^"'^, although including the ef- 
fect of relaxation brings it in line with the other volume- 
conserving constants. 



Solving the linear system of nine equations in the nine 
unknown force constants, we obtain the values listed in 
Table The volume force constant Cq is nearly 40% 
smaller than in pure Pt and 25% smaller than in a- 
Pt2Si. Nonetheless, the value is still sizeable and perhaps 
somewhat surprising given that we found no evidence of 
metallic-typ e bon ding in our analysis of the charge den- 
sity in Sec. [VP . The first- and second- neighbor Pt-Si 
radial force constants are quite large and nearly as large 
as the first-neighbor Si-Si force constant in pure Si. This 
result is consistent with the fact that we found only a 
small number of two- and three-center bonds for each 
Si atom in PtSi. This small number of bonds means 
that each bond is relatively strong, as is the case in pure 
Si, but in contrast to the situation in Q!-Pt2Si where the 
Pt-Si Ci force constant is more than a factor of three 
smaller. The fact that Ci and C2 in PtSi are still smaller 
than Ci is Si may be due to the fact that the bond angles 
in PtSi are considerable distorted away from the perfect 
tetrahedral angle. The Pt-Si C3 force constant in PtSi 
is approximately a factor of four smaller than Ci and C2 
which may be due in part to the correspondingly longer 
bond length. 

The Pt-Pt sixth-neighbor force constant Cg is larger 
than Ci in pure Pt which is likely due to the fact that 
this interaction contributes to the three-center bonds in 
PtSi. However, Cg is 30% smaller than the corresponding 
Pt-Pt C2 force constant in a-Pt2Si, reflecting the longer 
bond length and the lack of two-dimensional metallic 
sheets in PtSi. Although the seventh-neighbor Pt-Pt 
bond length in PtSi is only 0.03 A larger than the sixth- 
neighbor bond length, the seventh-neighbor bond does 
not participate in any three-center bonds and we found 
little evidence of any increase in the charge density. It is 
thus not surprising that C7 is more than a factor of four 
smaller than Cq. 

We find that the angular interactions are sizeable in 
PtSi, as they were in Q!-Pt2Si. However, in PtSi these in- 
teractions show a wider variation in magnitude, with if 23 
being more than an order of magnitude larger than K13. 
We can understand the variation in these Pt-Si-Pt force 
constants by looking at the sizes of the bond angles them- 
selves, if 13 corresponds to a bond angle of 71.09° which 
is very far from the perfect tetrahedral angle of 109.47°. 
The bond angle associated with if 22 is a lot closer, hav- 
ing a value of 94.64°, resulting in a larger force constant. 
The largest angular force constant is if 23 with the corre- 
sponding bond angle of 109.75° being nearly identical to 
the perfect tetrahedral angle. While the trend in the an- 
gular force constants in PtSi is understandable in terms 
of the deviation relative to the pure tetrahedral angle, 
the large magnitude of if 23 in comparison to the angu- 
lar force constants in pure Si and a-Pt2Si is unexpected. 
The angular interactions appear to be of greater impor- 
tance in PtSi than they were in a-Pt2Si. An attempt to 
fit the elastic constants of PtSi using a valence force field 
model including only radial interactions plus a volume 
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term resulted in nonsensical values for these force con- 
stants. A sensible fit was only achieved after including 
angular terms. 

We can now examine some of the trends in the elas- 
tic constants of PtSi using the calculated force constants. 
In particular, the Cauchy relations for an orthorhombic 
crystal are that ci2 — cee, C13 — C55, and C23 = C44.E3 
The expressions for the deviations from these Cauchy re- 
lations are as follows, 



(ci2 - cge) 



frozen 



-(Co - 0.1765 K22 - 0.0860 K23), 

V 



(35) 



(Cl3 - C55) 



frozen 



-(Co - 0.4475 Ki3 - 0.3453 X23), 



(C23 - C44 



, frozen 



(36) 



-{Cq - 0.3202 K22 - 0.4544 K23)- 

V 

(37) 



As in the case of a-Pt2Si the volume interaction makes a 
positive contribution to the deviations from the Cauchy 
relations while the angular interactions make a negative 
contribution. The geometry coefficients for the angular 
terms in Eqs. (35)-(|37|) are smaller than for Q;-Pt2Si in 
Eqs. (|3l| ) and (32), reflecting the smaller multiphcity of 
the bond angles in PtSi. This reduction is more than 
compensated by the larger magnitude of the force con- 
stants in PtSi, particularly K23. The volume per atom v 
is similar in the two silicides but the magnitude of Co is 
smaller in PtSi. The combined effect of the smaller Co 
and the larger ^^23 is that the deviations from the Cauchy 
relations in Eqs. (|35|)-(^7|) are still positive but approx- 
imately 30% smaller on average than in Q;-Pt2Si. This 
conclusion remains true for the relaxed elastic constants, 
although the specific numerical details are changed. For 
example, the larger geometry coefficients of ^22, and es- 
pecially ^^23, in Eq- ( p7| ) result in a very small deviation 
from the third Cauchy relation (C23 — 044)^™''®'^ for the 
frozen elastic constants. When relaxation is included C44 
drops by 29% while C23 increases by 8%, resulting in a 
significantly larger deviation. However, (ci3 — C55) be- 
comes much smaller so that on average the deviations 
are still approximately 30% smaller in PtSi. 

The requirements of mechanical stability in Eqs. (|^) 
and (M) constrain the elastic constants by requiring that 



C12 < 2(^11 



' C22), C13 < |(cii 



'C33), C23 < 2 ('^22 



So < 5(011 -I- C22 + C33), and Bq > i(ci2 -I- C13 



-C33), 

C23)- 

However, there are additional trends among the elas- 
tic constants. We have already noted that the volume- 
conserving elastic constants in PtSi are all smaller than 
those where the corresponding distortion does not con- 
serve volume. The predominant reason for this occur- 
rence is the presence of the positive volume contribution 
Co, just as it was in the case of pure Pt and in a-Pt2Si 



(see Fig. g). The positive deviations from the Cauchy re- 
lations in Eqs. (^5|)-(^ provide specific examples of this 
trend. We noted above that the relatively large value of 
Co in PtSi seemed surprising given the lack of evidence 
for metallic bonding in the charge density. In fact, it 
appeared that the chemical bonding in PtSi was much 
more similar to that in pure Si than in either a-Pt2Si or 
pure Pt. However, the trends in the elastic constants of 
PtSi, the positive deviations from the Cauchy relations 
and the smaller values of the volume-conserving elastic 
constants, are much more similar to those in the mate- 
rials that do exhibit direct evidence of metallic bonding, 
thus requiring a sizeable Co volume contribution in PtSi 
as well. This conclusion is not one that we would have 
reached based on the charge density alone, thus demon- 
strating the need for care when examining such quali- 
tative characteristics. By contrast, the analysis of the 
elastic constants using a valence force field model has 
allowed a more quantitative description of the chemical 
bonding. We note that the finding of both metallic and 
covalent components to the bonding in PtSi as well as 
Q!-Pt2Si indicates a strong similarity between these two 
materials and may also be connected with the fact that 
the heats of formation for the two are very nearly the 
same.Q 

Finally, we examine how the elastic constants of PtSi 
fit into the trends between the different materials studied 
here. The expression for the bulk modulus in PtSi is 



B, 



frozen 



1 



(Cll + C22 + C33 + 2ci2 + 2ci3 + 2C23)^''' 



i(co + ^Ci-Hic2 



1 



1 



+ T7;^3 + -^C6 + — C7 . (38) 
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Using the force constants listed in Table |^ we find that 
there is a roughly 50:50 split between the volume and 
radial contributions to Bq in Eq. (^8|) compared to an 
approximate 60:40 split in pure Pt and a-Pt2Si. We al- 
ready noted that the volume per atom v is similar in 
all three materials. Thus we see that the smaller value 
of Co is partially compensated by an increase in the ra- 
dial contribution, yielding a value of Bq which is only 
slightly smaller in PtSi than in a-Pt2Si, but still approx- 
imately a factor of two larger than in pure Si. We can now 
see that the nearly linear relationship between the bulk 
modulus and the atomic percent Pt, evident in Fig. |^, 
has a direct connection with the nature of the chemical 
bonding in these materials. Conversely, the fact that the 
volume-conserving elastic constants are similar in mag- 
nitude in all four materials demonstrates that they are 
less sensitive to the nature of the bonding. For example, 
C44 in Pt [Eq. (^0|)] arises purely from radial interactions 
while the split is 90:10, radial to angular, in Si [Eq. (pj|)]. 
In the two silicides [Eqs. (|2^) and (|33|)] the split is ap- 
proximately 60:40. Despite these variations in the split 
between radial and angular contributions and variations 
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in the individual force constants themselves, the volume- 
conserving elastic constants are all relatively small and 
similar in magnitude in all four materials. 



VI. SUMMARY 

We have carried out an extensive study of the chemical 
bonding and elasticity of two room-temperature stable 
platinum silicides, tetragonal a-Pt2Si and orthorhombic 
PtSi. We have investigated the trends in the calculated 
elastic constants, both the trends within a given material 
as well as between materials. The Cauchy relations, that 
ci2 = C66, ci3 = C55, and C23 = C44, apply to a crystal 
in which the interatomic interactions are purely radial. 
Real materials deviate from these relations and we find 
that in pure Pt as well as the two silicides the deviations 
are always positive (left hand side greater than right hand 
side) but in Si the deviation is negative. More generally, 
we find that in the metals the elastic constant expres- 
sions which correspond to volume-conserving strains are 
always smaller than those which correspond to strains 
which do not conserve volume. This also turns out to be 
true in Si with the exception of C12 which is less than C44 
(negative deviation from the Cauchy relation). However, 
the difference in magnitudes between volume-conserving 
and non- volume-conserving elastic constants is largest on 
average in Pt and gets smaller in the progression Pt — > 
a-Pt2Si —^ PtSi — > Si. In general, the volume-conserving 
elastic constants have similar magnitudes in all four ma- 
terials while the non- volume-conserving elastic constants 
follow this same progression. In particular, the bulk mod- 
ulus is found to be a very nearly linear function of the 
atomic percentage of Pt. 

We have analyzed the valence electronic charge density 
in order to gain insight into the nature of the chemical 
bonding in the silicides. In the case of Q;-Pt2Si we find 
striking evidence of a wide network of covalent three- 
center bonds, each involving a single Si atom and two Pt 
atoms. Each Si atom participates in 12 different three- 
center bonds. We also find evidence of two-dimensional 
metallic Pt (001) sheets which act to interconnect the 
network of three-center bonds. The Pt-Pt bond length 
in these two-dimensional sheets is very nearly the same 
as in pure fee Pt. The widely distributed nature of the 
bonding in Q;-Pt2Si appears to be closer in character to 
the pure metallic bonding in fee Pt than the covalent two- 
center bonds in Si. The trends in the elastic constants 
support this interpretation. PtSi also exhibits evidence 
of covalent Pt-Si-Pt three-center bonds in addition to 
more standard Pt-Si two-center bonds. Each Si atom 
participates in one three-center bond and two two-center 
bonds with the four Pt neighbors forming a very distorted 
tetrahedron. Two of the six corresponding bond angles 
are very nearly equal to the perfect tetrahedral angle but 
the other four angles vary from 71° to 132°. Qualitatively 
the bonding in PtSi appears much more similar to the co- 



valent bonding in pure Si than the metallic bonding in 
pure Pt but the trends in the elastic constants indicate 
that there are actually elements of both. The finding of 
strong Pt~Si covalent bonding in PtSi is-.consistent with 
the experimental study of Franco et alU in which they 
found spectroscopic evidence that the infiuence of the Pt 
6d orbitals extends throughout the entire valence band. 

We have constructed valence force field models for the 
two silicides as well as pure Pt and pure Si. These models 
provide a quantitative basis for understanding both the 
trends in the elastic constants and the various elements of 
the chemical bonding. We have included volume-, radial-, 
and angular-dependent contributions in the models. The 
volume-dependent contribution, which reflects the pres- 
ence of metallic bonding, turned out to be a crucial el- 
ement of the models. The presence or absence of this 
term and the magnitude of the volume force constant Co 
are predominantly responsible for the observed trend in 
the non-volume-conserving elastic constants as a func- 
tion of Pt concentration. In addition, the absence of this 
contribution in the volume-conserving elastic constants 
is largely responsible for the fact that these constants 
have similar magnitudes in all four materials. The vari- 
ation in the sign and magnitude of the deviations from 
the Cauchy relations is a specific example of these more 
general trends and is once again due primarily to the 
variation in the magnitude of Cq. The models also pro- 
vide explanations for differences in magnitude between 
specific elastic constants for a given material, such as the 
anomalously large value of cge in Q!-Pt2Si which we find 
to be closely connected to the three-center bonds in this 
material. 

In addition to providing explanations for the trends in 
the elastic constants, the magnitudes of the various force 
constants themselves provide a direct indication of the 
nature of the chemical bonding. The magnitude of the 
volume term provides an indication of the relative impor- 
tance of metallic bonding. This analysis demonstrated 
that there is an important element of metallic bonding in 
PtSi, despite the lack of direct evidence in the analysis of 
the charge density. This conclusion is required as a result 
of the specific values of the elastic constants in this mate- 
rial and would not have been possible based solely on the 
qualitative features of the charge density. Similarly, the 
magnitudes of the radial and angular force constants are 
directly connected to the importance of covalent bonds 
in the material. The trends in these constants confirm 
the general conclusions made on the basis of the charge 
density analysis. In addition, the conclusion that there 
are elements of both metallic and covalent bonding in a- 
Pt2Si as well as PtSi may be connected to the fact that 
the heats of formation for the two silicides are nearly 
the same. One general conclusion of this study is that 
the elastic constants contain a great deal of information 
about the nature of the chemical bonding in a material 
but since this information is not readily apparent, an 
analysis such as the one presented here is necessary in 
order to extract the information. We have attempted to 
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make the case here that an analysis in terms of valence 
force field models provides a convenient and fruitful way 
to analyze the elastic constants and their connection to 
the chemical bonding in a material. 

Our purpose in developing the valence force field mod- 
els described in this work was to provide a quantitative 
means for investigating the nature of the chemical bond- 
ing in the platinum silicides in comparison to pure Pt 
and pure Si and also to provide a more intuitive under- 
standing of the connection between the chemical bonding 
and mechanical properties of these materials. Nonethe- 
less, we can briefly consider the possibility that these 
models may be useful in carrying out future studies of 
silicide-silicon interfaces where first principles methods 
would be vastly more CPU-intensive. For example, de- 
pending on the growth conditions, the silicide thin film 
grown on a silicon substrate can be stabilized in an amor- 
phous phase. The only hope of treating such a structure 
would be to use a more efficient semi-empirical method 
such as a valence force field model. We believe that in 
general it should be possible to develop such a model 
given that our basic formulation includes the same fun- 
damental elements as in other successful models, such as 
the embedded-atom method and Tersoff potentials. One 
possible point of concern is the well-known fact that va- 
lence force field models in general tend to converge very 
slowly with respect to the number of interaction param- 
eters in the model. This issue would certainly need to be 
explored before any attempt was made to develop models 
that could be used in large-scale simulations. 
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Pt Si 




FIG. 1. Conventional unit cells of (a) body-centered tetrag- 
onal a-Pt2Si and (b) orthorhombic PtSi. The relevant lattice 
constant distances are illustrated in both cases. 
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FIG. 2. Trends in the elastic constants as a function 
of atomic percent Pt for pure cubic-diamond-phase Si, or- 
thorhombic PtSi, tetragonal a-Pt2Si, and fee Pt. The differ- 
ent curves correspond to the average values of different classes 
of the individual elastic constants, as specified in the legend. 
For example, in the case of the dotted-line curve labeled as 
Ci2, the line passes through |(ci2 4- C13 -I- C23) in the case of 
PtSi and through |(ci2 -I- 2ci3) for a-Pt2Si (cia = C23 for 
tetragonal crystals), while the open squares show the actual 
values of C12, C13 and C23, as appropriate for each material. 
The Co force constant curve is scaled by the inverse of the vol- 
ume per atom in order to be able to plot it on the same scale 
as the elastic constants. The significance of Co in connection 
to the elastic constants is discussed in Sec. 




FIG. 3. Superposition of free atom densities subtracted 
from the fully self-consistent crystal density for (a) cubic di- 
amond-phase Si and (b) fee Pt. In both plots there are 51 
contour levels plotted with pure black corresponding to —10 
and pure white to -f 10 millielectrons/bohr^ , as indicated in 
the scale bar. In (a) the x-axis is along [110] and the y-axis 
along [001] while in (b) the x-axis is [100] and the y-axis 
[010]. In both cases the calculations were carried out at the 
experimental equilibrium volume and only the density from 
the valence states was considered, excluding the density aris- 
ing from the core states. In the case of Si (a) the density was 
calculated at 76x101 grid points while for Pt (b) there were 
101x101 grid points. 
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FIG. 4. Superposition of free atom densities subtracted 
from the fully self-consistent crystal density for tetrago- 
nal a-Pt2Si. The same 51 contour levels and gray-scale 
are used as in Fig. ^. The two-dimensional Pt-Pt sec- 
ond-nearest-neighbor metallic sheets are shown in (a) with 
the X-axis along [100] and the y-axis along [010]. In (b) we 
show the three-center Pt-Si-Pt covalent bonds involving two 
second-neighbor Pt atoms, with the x-axis along [llO] and the 
y-axis along [111]. The second set of three-center Pt-Si-Pt 
covalent bonds involving two third-neighbor Pt atoms is illus- 
trated in (c) with the x-axis along [100] and the y-axis along 
[001]. All three calculations were carried out for the exper- 
imental equilibrium structure and only the density from the 
valence states was considered. The density was calculated at 
101x101 grid points in (a), 151x201 grid points in (b), and 
101x151 grid points in (c). 




FIG. 5. Superposition of free atom densities subtracted 
from the fully self-consistent crystal density for orthorhom- 
bic PtSi. The same 51 contour levels and gray-scale are used 
as in Fig. |^. In (a) we show the three-center Pt-Si-Pt covalent 
bonds with the x-axis along [lOO] and the y-axis along [001]. 
The covalent Pt-Si bonds which connect atoms in adjacent 
b-axis planes are shown in (b) with the x-axis approximately 
along [304] and the y-axis along [010]. Both calculations 
were carried out for the experimental equilibrium structure 
and only the density from the valence states was considered. 
The density was calculated at 101x101 grid points in (a) and 
201x65 grid points in (b). The leftmost Pt atom in (b) with 
the label contained inside a dotted circle is not actually lo- 
cated in the plane of the plot but is close enough that its 
influence can still be seen. 



TABLE I. Equilibrium theoretical (from Ref. |6|) and experimental lattice constants (in a.u.) and internal structural param- 
eters (for PtSi). 
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expt. 
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a-PtaSi 
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VSi 


Ref. 
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0.9977 


0.1919 


0.1782 


0.5841 


i 




expt. 


0.9956 


0.1922 
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0.583 





TABLE II. Elastic constants of Pt and Si. The flrst principles calculations, described in Ref. |6|, were carried out at the 
theoretical self-consistent lattice constants of apt — 7.403 a.u. and agi — 10.22 a.u. The theoretical value of C44 in parentheses 
for Si is the "frozen" value obtained without allowing for internal relaxation. The bulk modulus is calculated from the elastic 
constants as Bo = |(cii + 2ci2). Experimental values are extrapolated to K. All values are in units of GPa. 



Pt TheoryE 


Pt Expt. 


31 


Si Theory 


i Si Expt 


16 


Cll 


346.8±0.5 


358 




163.45±0.03 




165 




C12 


262.7±0.3 


254 




62.13±0.02 




63 




C44 


87.5±0.3 


77 




79.85±0.02 (108.6) 


79.1 




Bo 


290.8±0.3 


288.4 




95.90±0.02 




97.0 
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TABLE III. First principles elastic constants of a-Pt2Si and PtSi from Ref. |6[ Calculations were performed at the theoretical 
self-consistent lattice constants (Table 0). "Frozen" refers to keeping the atoms held fixed at the positions determined solely from 
the strain tensor and, in the case of PtSi, with the internal structural parameters held fixed at their theoretical self-consistent 
values (Table |). "Relaxed" indicates that a relaxation of the atomic positions was carried out, including a relaxation of the 
PtSi internal structural parameters. Parentheses in the case of the relaxed a-Pt2Si elastic constants denote values where no 
internal relaxation was necessary because of symmetry constraints (small differences with the frozen values come from using 
a slightly more stringent convergence criterion on the energy). The bulk modulus is calculated from the elastic constants as 
Bo = |(2cii + C33 + 2ci2 -I- 4ci3) for Q:-Pt2Si and Bo = ^(cii -I- C22 + C33 + 2ci2 + 2ci3 + 2c23) for PtSi. No experimental data 
is available for either material. All values are in units of GPa. 
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327.5±1.2 


298.2±1.2 


C22 






313.8±0.0 


269.3±0.8 


C33 


297.5±0.5 


(298.0±0.4) 


345.9±0.1 


308.0±0.6 


C12 


225. Oil. 2 


239.6±1.0 


157.7±0.6 


156.4±0.8 


Cl3 


169.3±0.9 


(169.4±0.8) 


162.9±0.6 


132.2±0.7 


C23 






153.4±0.1 


165.1±0.6 


C44 


75.4±0.3 


62.7±0.5 


141.3±0.3 


100.1±0.4 


C55 






113.1±0.1 


104.5±0.1 


C66 


169.5±5.2 


(169.3±5.2) 


74.2±0.2 


66.3±0.4 


Bo 


235.4±0.6 


(235.5±0.5) 


215.0±0.2 


198.1±0.3 



TABLE IV. Force constants of valence force field models [Eq. (|8|)] for fee Pt, tetragonal a-Pt2Si, orthorhombic PtSi, and 
cubic-diamond phase Si. The Cq force constant represents the volume-dependent interaction, each of the remaining d is a 
radial force constant for the ith nearest-neighbor bond, and Kij is an angular force constant for the bond angle between the 
ith and jth nearest-neighbor bonds. All of the force constants are in units of eV. 





Pt 


a-PtaSi 


PtSi 


Co 


16.54 


13.88 


10.38 


Ci 


16.10 


13.58 


42.58 


C2 




26.03 


48.26 


C3 




8.39 


10.90 


Ce 






18.32 


C7 






3.91 


A"ii 




5.06 








1.87 










1.29 


K22 






7.62 


K23 






15.01 
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